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SUMMARY 


The valid and invalid application of the Mikhailov criterion to linear, time -invariant 
systems with time delays is discussed. The Mikhailov criterion is a graphical procedure 
which was initially developed to examine the stability of linear, time-invariant systems 
with no time delays. For these systems, there are two equivalent formulations of the 
criterion. Attempting to apply the second formulation when there are time delays in the 
systems can lead to erroneous results, as shown by an example. However, the first formu- 
lation remains valid for time-delayed systems of the retarded type, with the understanding 
that the Mikhailov curve need not necessarily always rotate in the counterclockwise direc- 
tion for a stable system. At present, application of the Mikhailov criterion to neutral sys- 
tems has not been justified. 


INTRODUCTION 

The Mikhailov criterion is used to examine the stability of dynamical systems that 
are described by linear ordinary differential equations with constant coefficients. 
Kashiwagi (ref. 1) indicates that this criterion can be used with only limited success when 
there are constant time delays in the system. Yet recently, Chen and Tsay (ref. 2) have 
shown that a modified Nyquist criterion, which is similar to the Mikhailov criterion, does 
have valid application to time-delayed systems. Why is there a discrepancy between the 
earlier and recent works? The purpose of this paper is to remove the confusion related 
to the valid application of the Mikhailov criterion to time-delayed systems. 

SYMBOLS 


A,Aj,Bj real constants 

i imaginary unit, 

j integer index 


k 


real number 



L{s) characteristic polynomial for system with no time delays, or characteristic 

quasi-polynomial for systems with time delays 

LiIsjjLgCs) parts of L(s) shown in equations (42) and (43) 

M number of roots of L(s) with positive real parts 

N order of system with zero time delays 

s complex variable 

sj jth root of L(s) = 0 

X(w) real part of L(icc) 

Y(t o) imaginary part of L(ico) 

y angle in figure 2 

£ damping parameter 

0(a>),0,(a>) change in arguments of L(iw) and ico - s-, respectively, as w increases 

J J 

from zero 

8,6- total change in arguments of L(iu>) and ico - s-, respectively, as w varies 

J J 

increasingly from zero to infinity 
o real part of s 

asymptote of real part of large modulus roots of L(s) 
t, Tj constant real-time delays 

cu imaginary part of s 

ANALYSIS 

By referring to the complex geometry involved, Popov (ref. 3) gives insight into the 
Mikhailov criterion for a system with no time delays. Kashiwagi used Popov as a basic 
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reference. Hence, in order to appreciate the earlier thoughts on the Mikhailov criterion, 
it is necessary to first examine the criterion for no time delays as viewed by Popov. 

Mikhailov Criterion (No Delays) 

Let the characteristic polynomial of a linear system be expressed as 

L(s) = AqS^ + A^s^ + • • • + Aj s j_^s + A-^ ^Aq > 0) (1) 

with real coefficients, or in factored form as 

L(s) = A q (s - s^ (s - s 2 ) . . . (s - s N ) (2) 

where the complex roots appear in complex conjugate pairs. With s = icu, equation (1) 
can be written as 


L(iw) = X(w) + iY(cu) 

where 

X(w) = A n - A n _ 2 w 2 + A n _ 4 o> 4 - . . . 
Y(w) = Aj^jW - A n _ 3 cu 3 + . . . 


(3) 


(4) 


and equation (2) becomes 

L(iw) = A 0 ^iw - sj fiu - s 2 ) . . . (icu - s N ) (5) 

Now consider the plot of L(iw) in the complex plane as w is increased from 
0 to oo, as depicted hypothetically in figure 1. This is the so-called Mikhailov curve. 

The argument of L(iw) is simply the sum of the arguments of the linear factors of 
L(ico) in equation (5), that is, 

N 

(iw - sA (6) 

j=l V ' 


L(iw) = y 
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Let 0(co) and 0-(co) denote the changes in the arguments of L(ico) and ico - s-, 
respectively, as co increases from zero. Then, 

N 

0(w) = ^ 0-(w) (7) 

j=l 

(The argument of the real number A Q > 0 is zero.) Furthermore, if 9 and 0. denote 
the total change in the arguments of L(ico) and ico - Sj, respectively, as co varies 
from 0 to °°, then 

N 

~ e =l~° j (8) 

j=l 

Consider the following two special cases to gain insight into the contributions of 9- 
to 9 in equation (8). 

Negative real root .- The change in the argument of the factor ico - s^, when Sj is 
a negative real number, is examined by using figure 2(a). It is geometrically clear that 
as co varies from zero to infinity 

-- lim 0 1 (co) = J (9) 

co— 00 z 

Complex conjugate roots .- Consider the angular change in the arguments of the two 
factors ^ico - s 2 ) and ^ico - s^, when s 2 and s^ are complex conjugate roots with 
negative real parts, as indicated in figure 2(b). Geometrically, it can be seen that for 
co = 0, the argument of the complex vector ico - s 2 is -y. As co increases from zero 
to infinity, the change in the argument is 

0 9 = lim 0 9 ( co) = J + y 
z co— °° z ' 2 


Similarly, for the vector ico - S3, 


(10) 



Hence, 


e 2 + e 3 = I + f 


( 12 ) 


It is indicated by equations (9) and (12) that each root of L(s) with a negative real 
part contributes ^ to the argument of L(iw) in equation (8). By similar geometry, it 
can be seen that any root of L(s) with a positive real part will result in an argument 
change in L(iw) of as co goes from 0 to °°. Popov (ref. 3) also discusses the 
effects of a zero root, purely imaginary roots, and an infinite root, but these cases are 
not presented here. 

Suppose there are M roots of L(s) in equation (1) with positive real parts, and 
N - M roots with negative real parts; then it is not difficult to infer geometrically that 


e = (N-M)|-M|=(N-2M)| (13) 

In order for the systems under consideration to be stable, all roots of L(s) must 
have negative real parts. In this respect, Popov (ref. 3) writes: "For stability of an 
nth-order linear system it is necessary and sufficient that the Mikhailov curve plotted for 
the characteristic equation of the given system pass through n quadrants in succession 
counterclockwise, circling the origin of coordinates." Thus, since n = N in this paper, 
L(icu) completes a rotation by the angle 


e = n| (14) 

which is equation (13) with M = 0 (no roots with positive real parts). This is the first 
formulation of the Mikhailov criterion. 

The roots Sp s 2 , and Sg of L(s) are shown in figure 2 and have negative real 
parts. The associated vectors iw - Sp iw - s 2 , and iw - always rotate counter- 
clockwise as u) increases; or equivalently, 6 p 0 2 , and 6 3 always increase. There- 
fore, from equation (7), if all roots of L(s) have negative real parts, it is not possible 
for the vector L(iw) to ever rotate by any amount in the clockwise direction; that is, the 
vector L(ico) continually rotates in the counterclockwise direction as w increases 
from 0 to °°, and 0(u>) approaches N| (eq. (14)). This leads to the second formulation 

of the Mikhailov criterion, which states that a necessary and sufficient condition for a 
linear, time -invariant system with no time delays to be stable (characteristic roots with 
negative real parts) is that the real part, X(w), and the imaginary part, Y(w), of L(ico) 
in equation (3) alternately vanish a finite number of times. 
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Transcendental Characteristic Equation 

Previous works (refs. 3 and 4, for example) interpret the Mikhailov criterion in 
terms of a polynomial characteristic equation. Difficulty occurs in trying to apply this 
interpretation to systems with time delays (ref. 1). 

Suppose the system has time delays so that the characteristic function is, for 
example, the transcendental characteristic quasi-polynomial 

L(s) = s 2 + 2£se -TS + 1 (15) 

where r is a constant time delay. Unlike equation (1) which has exactly N roots, 
equation (15) has an infinite number of roots because of the exponential term introduced 
by the time delay. 

For illustration, let r = 6.2 and £ = 0.2 in equation (15). In this specific case, 
all the roots of equation (15) have negative real parts, making the system stable (ref. 1). 

The Mikhailov curve corresponding to equation (15) is shown in figure 3. Notice 
that the curve moves in both the clockwise and counterclockwise directions even though all 
the roots have negative real parts. This could never happen if L(s) were a polynominal 
with negative real parts. To approximate the curve in figure 3 very closely with a poly- 
nomial would require that the approximating polynomial have an unstable root (positive 
real part) in order that the approximating curve ever move in a clockwise direction over 
any portion of the solution. 

Clearly, from figure 3, X(w) and Y(o>) do not alternately vanish; and, conse- 
quently (by counterexample), the second formulation of the Mikhailov criterion definitely 
does not apply in general. However, as noted in the section entitled "Mikhailov Criterion 
for Delayed Systems," with proper interpretation, the first formulation of the Mikhailov 
criterion can be applied. 


Mikhailov Criterion for Delayed Systems 

Recently, Chen and Tsay (ref. 2) derived a modified Nyquist stability criterion which, 
in the context of this paper, is similar to the first formulation of the Mikhailov criterion. 
The criterion was shown to hold for characteristic functions L(s) for which: 

(1) L(s) has no purely imaginary roots. 

(2) L*(s) = L(s*), where * means complex conjugate. 

(3) L(s) acts like As“k as s — °° in the right half of the s-plane; that is, 
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(a g 0) (16) 


lim = A ^ 0 

S — °o S - ^ 

where A is a nonzero constant and k is an integer. 

Chen and Tsay (ref. 2) show that, as o> increases from 0 to °°, the vector L(io>) 
rotates by the amount 

?=(-k-2M)| (17) 

where M is the number of roots with positive real parts. 

Equation (17) indicates a stable system if and only if M = 0; that is 

e = -k| (is) 

It is important to note that there is no requirement that 6 result only from counterclock- 
wise rotations. 

With respect to the characteristic polynomial in equation (1), k = -N in equa- 
tion (16), and equation (18) becomes 

o = n| (19) 

which is the same as equation (14). Now, from Popov's view of the Mikhailov criterion, 
it is clear that 6 in equation (19) will, indeed, consist only of counterclockwise rotations. 
However, this is not explicit in the derivation (ref. 2) of equation (18). 

Retarded systems .- Consider a system with a characteristic quasi -polynomial of 
the form 


N N-l _ Tg 

L(s) = ^ ^ BjS^e (20) 

j=° 3=o 

where A^ ^ 0 and Tj g 0 (ref. 5). The dominant term in equation (20) as s — °° is 
N 

Ajq-s . To show this, consider the magnitude of the ratio of any remaining term in equa- 
tion (20) relative to A^s 1 ^. Thus, note that for j ^ N, 
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( 21 ) 



v 

lim 

l s b°° 

A n s N 

since N - j > 0. 

Also, 


B-S'i 

S 

lim 

| S | — °o 

A S N 

a n s 


lim 

l s h°° 


h. 

a N| 


lim 

\s\~°o 


s 


I N— j 


= 0 


B. 




~ T f 

sl N "j 


= 0 


( 22 ) 


where o- is the real part of s. Hence, 

lim = A n 4 0 (a £ 0) (23) 

s — oo s iN 

for the retarded system, and it follows from equation (16) that k = -N. Therefore, from 
equation (18), the system is stable if and only if 


? = N — (24) 

2 

For example, suppose the retarded system has the characteristic quasi-polynomial in 
equation (15). Then, equation (24) indicates that L(ico), which is shown in figure 3, should 
rotate by the angle 2 = n for stability. The real part of L(ico) in figure 3 is 

increasing much faster than its imaginary part; therefore, the change in the argument of 
L(ico) as <o varies from 0 to °° appears to be approaching 7 r. This conjecture is veri- 
fied algebraically as follows. With s = ico, equation (15) becomes 

L(iu>) = (l - o>2 + 2£co sin cor) + (2£co cos wr)i (25) 

The real part of L(ico) in equation (25) is 

X(co) = 1 - co 2 + 2£co sin cut (26) 

Notice in equation (26) that if 

co 2 > 1 + 2£co sin cor (27) 
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then X(co) <0. Certainly, equation (27) will be true if 

co 2 > 1 + 2£co (28) 

The right-hand side of equation (28) dominates the left side as co approaches 0. How- 
ever, as co increases, the left side will eventually dominate the right, and equation (28) 
will be valid. Equation (28) is valid for all co > co m , where co m is the largest positive 
real root of the polynomial equation 


co 2 - 2£co -1 = 0 ( 29 ) 

and is given by 

"n, = ? + V? 2 + 1 (30) 

For the specific example being considered, £ = 0.2 and t = 6.2. Thus, equations (26) 
and (30) become, respectively, 


X(co) = 1 - co 2 + 0.4co sin 6.2co 


(31) 


co 


m 


1.2 


(32) 


The real and imaginary parts of L(ico) are parametric equations in co, with co 
increasing along the Mikhailov curve in figure 3. After co increases to co = co m = 1.2, 
the corresponding X(co) value is X(co m ) = 0.056. Thereafter, co > co m and corre- 
sponding values of X(co) will remain negative. In other words, if continued, the 
Mikhailov curve in figure 3 will remain in the left half plane (second and third quadrants). 
Since the curve in figure 3 does not encircle the origin, the change in the argument of 
L(ico) as co varies from 0 to °° lies in the interval 


J < e < — 
2 2 


(33) 


Since the argument of L(ico) in figure 3 is zero at co = 0, 6 ( co) = arg (L(ico)]. 
Hence, 


cos 0(co) = 


1 - co 2 + 2^co sin cor 


cos c ot)^ + (l - co^ + 2^co sin cor) 


(34) 
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from which it follows that 


lim cos 0(co) = -1 

CO — 00 

Equation (35), along with inequality (33), implies that 0(w) - it as co - «>, so that the 
system is stable for t= 6.2 and £ = 0.2. 

Kashiwagi (ref. 1) was dealing with retarded systems when he stated that the 
Mikhailov criterion could be used with only limited success when there are constant time 
delays in the system. He draws this conclusion from the second formulation of the 
Mikhailov criterion, which is associated with characteristic polynomials. Indeed, as 
substantiated by the example in this paper, the second formulation does not hold in general 
for retarded systems; however, the first formulation remains valid if it is recognized 
that the Mikhailov curve need not always rotate counterclockwise. The only requirement 
is that the curve complete a rotation of N| (eq. (24)). This appears to be the basic 
reason for the discrepancy between earlier and recent works. 

The first formulation of the Mikhailov criterion can be used to determine if a 
retarded system is stable for a fixed set of time delays. In programing the criterion on a 
digital computer, it is not necessary to express the characteristic equation in simplified 
form, for example, like equations (1) or (3). 

The Mikhailov criterion cannot be used readily to generate stability boundaries. 

In this respect, however, parameter-plane methods (refs. 6 and 7, for example) are per- 
haps more appropriate. 

Neutral systems.- A system having the highest ordered derivative to appear both 
with and without a time delay is called a neutral system. A class of neutral systems has 
the following characteristic quasi-polynomial (ref. 5): 

N . N . _ T . S 

L(s) = ) A.s ] + 2 B. s ] e 3 ( 36 ) 

fo j=0 

where A N 4 0, B N 4 0, r. g 0 for j ^ N, and t n > 0. (For A N 4 0 and B N = 0, 

the neutral system reduces to a retarded system, provided there is at least one Tj 4 0 
for j 4 N.) These systems are stable if all the roots of equation (36) have negative real 
parts and are not asymptotic to the imaginary axis (ref. 8). The stability is uncertain if the 
roots are asymptotic to the imaginary axis. However, it appears that a system with such 
roots would be too lightly damped. In addition, because of the uncertainty in system coef- 
ficients, it seems desirable not to have a root too close to the imaginary axis. Within the 
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context of this paper, this condition is simply considered unsatisfactory. The system is 
unstable if any root has a positive real part. (See ref. 9.) A practical example of a 
neutral system is examined in reference 10. 

The stability of neutral systems has been determined recently by the extended 
T-decomposition method (refs. 11 and 12). In this section, the possibility of also applying 
the first formulation of the Mikhailov criterion to the class of neutral systems is 
examined. 

The characteristic equation for the neutral system is L(s) = 0, which can be 
expressed as 


N-l / 


A N sN = B N sNe 




+ 


A i + B i e 


-T jS 


j=0 


All roots of equation (37) must satisfy the relation 


(37) 



(38) 


which is obtained by applying properties of the absolute value to equation (37). Implicit in 
inequality (38) is the requirement that as |s| — 00 (he., the modulus of the roots increases 
without bound), 

i *}”„ (m - N e Tif ) = 0 <39) 

From equation (39), a becomes arbitrarily close to 


cr 


OO 




(40) 


That is, the real parts of the characteristic roots form a sequence which approaches cr^. 
As long as ^ 0, the roots will not be asymptotic to the imaginary axis. However, it 
still remains to be determined whether there are any roots with positive real parts 
(unstable) . 
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Three distinct cases are of interest, namely: 


Case 1: — — <1 or cr >0 

p °o 

b n 


(unstable) 


Case 2: — — =1 or a =0 

TD OO 

B N 


(unsatisfactory) 


N 

Case 3: — - >1 or a <0 

■R 00 

b n 


(either stable or unstable) 


In equation (40), > 0. Thus, for Case 1, cr^ > 0 and the neutral system is unstable. 

For Case 2, = 0, the roots are located too close to the imaginary axis, and the system 

is considered unsatisfactory. Finally, in Case 3, < 0 and the neutral system is 

either stable or unstable, being stable if and only if all roots have negative real parts. 

The remainder of this paper is concerned with Case 3 of the neutral system. If it 
could be shown that L(s) in equation (36) satisfies equation (16), then the previous 
results of Chen and Tsay (ref. 2) could be used to justify the application of the Mikhailov 
criterion. However, as the following analysis shows, this justification was not possible. 

For convenience of discussion, define 


L(s) = L x (s) + L 2 (s) 


(41) 


where 


L x (s) 



+ B N e 



N 


and 


(42) 


L 2 (s) = 



+ 


b j 6 



(43) 


It will be shown that L^(s) is the dominant term of L(s) as s — °° in the right 
half of the s-plane (a g 0). Inequality relationships involving the magnitudes of L 1 (s) 
and L 2 (s) are obtained from equations (42) and (43) as 
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(44) 


L 


1 



> 




e 


- T N a 


N 


and 



(45) 


If the coefficient of js|^ in equation (44) is bounded away from zero (greater than some 
positive number) and the coefficient of each |sp in equation (45) is less than some finite 

Lj(s 

half of the s-plane. There is no conflict with equation (39) because all but finitely many 
roots of the characteristic equation are in the left half of the s-plane for Case 3. 


positive number, then eventually 


)| will dominate L 2 (s)| as s — °° in the right 


For 0 < a < oo and since 
bounded away from zero as 




the coefficient on the 



term is 


0 < I a n - B N = 


a n ' b n 


- t N CT 


(46) 


so it is greater than the positive finite number |A^| - |B^|. 

For 0 < a < the coefficient of the |sp term is bounded above as 



+ 





+ 



(47) 


so it is less than the finite number 
is L^(s) as s — °o with ct 5 0, or 


A- 

+ 

B- 

J 


J 


Hence, the dominant term of 


L(s) 


lim 

s— 


L(s) 


s ^ A N + B N e 


■ t n s 


= 1 


(a ^ 0) (48) 


Because of the exponential term in the denominator, equation (48) is not in the form of 
equation (16). Hence, at this point, application of the Mikhailov criterion cannot be justi- 
fied for neutral systems on the basis of the results in reference 2. 
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CONCLUDING REMARKS 


The valid and invalid application of the first and second formulations of the Mikhailov 
criterion to linear, time -invariant systems with time delays has been examined. It has 
been demonstrated by an example that the second formulation of the Mikhailov criterion 
definitely should not be applied to systems with time delays. The first formulation, 
however, is directly applicable to delayed systems of the retarded type. 

These facts appear to be the basis of the discrepancy in earlier and recent works. 
The statement that the Mikhailov criterion had limited validity to time-delayed systems 
in the earlier work was based only on the application of the second formulation to retarded 
systems whereas recent work is equivalent to the first formulation, which does have valid 
application for such systems. 

At present, application of the Mikhailov criterion to neutral systems has not been 
justified. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 22, 1978 
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